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Abstract 


The origin of continents is not well understood, perhaps because each of the two usual ways of address- 
ing the problem has methodological difficulties. Isotope geochronology has enabled the condition of the 
early Earth to be reconstructed from the direct evidence of successively older rocks, working backwards 
in time. But, since a point is reached where there are no older rocks available, much speculation about the 
immediately pre-geological history of the Earth (the period when its identity was well established but of 
which no direct trace remains in the form of existing surface rocks) has come from planetology; and its 
approach has been, properly, to work forwards in time, starting from some preferred model of accretion. 
This paper is an attempt to bridge these two approaches, by focusing on what characteristics of the 
earliest rocks actually need to be explained by the planetological models. 


The restriction of these oldest rocks to a small number of core areas within early cratons, of which the 
Kaapvaal Craton (KC) of southern Africa and the Pilbara Craton (PC) of Western Australia are taken as 
typical, is first noted. The closely parallel evolution of both of these Cratons starts at about the same time, 
with the sudden appearance of vast quantities of evolved silicic magma (now orthogneisses), accompa- 
nied by the deposition of thick successions of mixed, mainly volcanic, rocks differing in only a few 
respects from those that have continued to come from the Earth’s interior; although the early Earth must 
have been hotter, this is not reflected in the metamorphism of the oldest strata. Following this abrupt 
beginning, these protocontinents grew steadily through continuous, but diminishing, magmatic activity, 
manifest as both silicic plutonism and bimodal eruption. 


All of these features are consistent with the following simple model of the early Earth. A completely 
molten immediately pre-geological Earth had, outside its metallic core, a pattern of huge, polygonal, 
convection cells with ascending central plumes and descending peripheries; there were possibly as few as 
12 such cells, forming a pattern of irregular pentagons, with 20 main centres of convective descent 
(COCDs) at their triple junctions. Before the continents formed, the thin skin of cooled and solidified 
material, formed on the surface of each cell during radial surface flow, was carried down at the margins 
and remelted as it descended. Such melting produced, by well known petrogenetic processes, a variety of 
differentiated magmas which during rapid early convection were carried down and remixed with their 
parent materials. As the whole Earth cooled, and the rate of convection slowed, a point was reached 
where the buoyancy of these silicic magmas enabled them to rise faster than the descending columns 
below the COCDs, resulting in the appearance at the surface of a plug of orthogneiss. 


The subsequent geological development of the KC and PC, jointly taken as archetypes of the early 
continents, fits comfortably with a concept that later continental growth was a simple response to mainte- 
nance of the same large-scale convective system which remained in force long enough to establish, below 
the oldest parts of the continents, the very deep roots for which there is increasing evidence. At some 
stage of lateral growth, possibly associated with both a change of global convective patterns and the 
development of a layered mantle, the early continents became unstable and were subject to break-up and 
relative displacement over the Earth’s surface. The existing areas of oldest rocks are not simply the few 
surviving remnants of some early (but not earliest) continental crust that covered the whole Earth, but 
really do represent both the initial and sole centres of continental growth. 
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Abstract 


The magnitude of the mass of an electron neutrino is still uncertain, yet this parameter is of vital 
importance to our understanding of cosmology, particle physics and neutrino oscillations. The 
process of double beta decay offers a means of testing grand unification theories as it may yield 
limits for the electron neutrino mass and indicate the degree of lepton conservation in elementary 
particle reactions. A number of measurements of the half-life of double beta decay of Te) 
1012%Xe and “Se/*Kr have been obtained, but the values are almost certainly inaccurate because of 
leakage of the gaseous daughters over geological time. This may be overcome by measuring the 
double beta decay half-lives of nuclides which have solid daughter products. Improvements in 
solid source mass spectrometry now allow such measurements to be successfully attempted. The 
results of one such experiment will be described and other double beta decay systems which may 


yield experimental results suggested. 


Introduction 


On December 4' 1930, the Austrian physicist 
Wolfgang Pauli suggested that a particle might be re- 
sponsible for carrying off the varying amounts of energy 
that seemed to be missing in beta decay. However, over 
two decades were to elapse before neutrinos were finally 
discovered by Reines & Cowan (1956). It has since been 
shown that there are different types of neutrinos - elec- 
tron neutrinos which are emitted during beta decay, 
muon neutrinos which are associated with charged pion 
decay, and tau neutrinos. 


The only way that neutrinos can interact with normal 
matter is through an extremely feeble force known as the 
weak interaction. Thus, in the “standard” model of par- 
ticle physics, neutrinos are mass-less particles that do not 
disintegrate. There is some experimental evidence (albeit 
controversial) that neutrinos can oscillate from one form 
to another (Athanassopoulos et al. 1995). This can only 
occur if these neutrinos possess a mass. In 1987, super- 
nova SN 1987A exploded in the Large Magellanic Cloud. 
From the spread in the arrival times of the burst of neu- 
trinos emitted by SN 1987A, it could be calculated that 
the mass of the electron neutrino is less than 25eV. 


One of the most successful developments of modern 
science has been the rapid evolution of nuclear astro- 
physics to what is now a mature scientific study. A 
study of nuclear astrophysics has enabled us to unravel 
many of the secrets of nucleosynthesis ~ those nuclear 
constraints which determine the abundances of the 
chemical elements. The isotope abundances of the ele- 
ments, which we measure so precisely by mass spectrom- 
etry today, have been fashioned in the stars and interstel- 
lar medium over billions of years. Some are vestiges of 
the Big Bang itself. 
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The three basic process involved in synthesising the 
heavy elements (Z>31) are the slow (s) and rapid (r) neu- 
tron capture processes, and the proton (p) process, a re- 
action that synthesises the neutron-deficient isotopes. 
The s-process synthesises nuclides along the valley of 
nuclear stability during the red giant phase of a star's 
evolution, whilst the r and p processes occur on a much 
shorter time scale in cataclysmic, non-equilibrium envi- 
ronments, such as supernova. De Laeter (1990) has de- 
scribed these nucleosynthetic processes in greater detail. 


Figure 1 shows a section of the Chart of the Nuclides 
in the vicinity of tin. The r-process nuclides mimic fis- 
sion product decay chains and produce neutron-rich iso- 
topes, in some cases far from the valley of stability. From 
the nuclear stability of these r-process even-even nu- 
clides, it can be shown that they cannot be absolutely 
stable, yet the fact remains that they have sustained their 
existence in measurable quantities over geological time. 
There is therefore a need for a nuclear decay mechanism 
to exist that is sufficiently slow to enable these nuclides 
to survive in significant quantities over a period of many 
billion years. 


Double beta decay is just such a mechanism in that its 
half-life is in excess of 10” years (Mayer 1935). The 
nuclear systematics of double beta decay is applicable to 
even-even nuclides on the neutron-rich side of the valley 
of nuclear stability, thus enabling survival of these nu- 
clides over geological time scales. 


Double beta decay 


Double beta decay is the rarest nuclear phenomenon 
in Nature. In double beta decay, a nucleus undergoes a 
transmutation from one element to another such that two 
electrons are emitted rather than the single electron emit- 
ted, in the more commonly observed single beta decay. 
All of the parent and daughter isotopes involved in 
double beta decay are even-even nuclei. The pairing 
force acting between like nucleons is responsible for the 
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Figure 1. A section of the Chart of the Nuclides in the vicinity of tin. The arrows represent nuclides produced by the r-process, which 


indicates the various modes of nucleosynthesis. 


increase in the binding energy of these nuclei relative to 
the odd-odd isotopes, thus preventing a single beta decay. 


The long standing question in double beta decay is 
whether the electron neutrino should be described by a 
Dirac or by a Majorana field (Primakoff & Rosen 1981). 
In the Dirac formalism, the neutrino has a distinct anti- 
particle such that: 


(A, Z) 3 (A,Z+2) + 2e + 20, 


However, if the neutrino is a Majorana particle then 
the neutrino and anti-neutrino are indistinguishable, so 
that double beta decay can occur with the net emission of 
no neutrinos such that; 


ON Fi ON eR OY 


Grand Unified Theories of particle physics allow the 
neutrino to have a finite mass, thus allowing the possibil- 
ity of neutrino-less double beta decay. However, double 
beta decay systems have such long half lives (10% - 10% 
years; Doi et al. 1985), that enormous experimental con- 
straints are imposed on the detection of double beta de- 
cay daughter products. 


There are two approaches to measuring the half-life of 
a double beta decay reaction. In the direct detection 
method, coincidence counting is used to measure the ra- 
dioactivity at the energy value predicted from precise 
atomic mass measurements (Dyck et al. 1990). The disad- 
vantage of the radioactive counting method is that intrin- 
sically faster radioactivities, such as those associated with 
uranium or thorium decay chains, produce significant 
backgrounds which may mask the beta decay product. 
Instrumental experiments are designed to detect the elec- 
trons as they are emitted and are thus capable of not 
only measuring the lifetime of the double beta decay, but 
also the electron energy spectrum. Doi et al. (1988) have 
summarised a number of experimental results for the 
half-life of double beta decay un Ge, “Se, ""Mo, '*Xe 
and “Nd conducted by a number of research groups 
(Table 1). The half lives range from >6 x 10' years to 
>1.4 x 10” years. 


One of these direct detection experiments was con- 
ducted by the University of California, Santa Barbara and 
Lawrence Berkeley Livermore, 600m underground with 
eight 0.9kg Ge detectors. The detectors are inside a cav- 
ity formed by 10 blocks of 15 cm thick sodium iodide 
crystals providing an active volume with a 30 keV 
threshold. The sodium iodide crystals are in turn inside 
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Table 1 


The data of the half-lives of the neutrinoless double-beta decay 
with Majorana emission for various nuclei. (Doi 1988 and pri- 
mary references therein). 


Experimental group Half-life (yr) 


“Ge Osaka 22102 
Moscow (ITEP) 20105 
Batelle-South Carolina 6 (+1) 10” 
Caltech-SIN-Neuchatel 120102 
Santa Barbara-LBL >14 10% 
8Se Irvine >4.4 10% 
1WMo Osaka >6 10" 
Irvine >75 10" 
136Xe Milano ECHO 
Moscow (INR) S10 107 
MON d Moscow (INR) >7.0 10" 


a pure Pb shield of 20 cm thickness. The detectors are 
kept at liquid-nitrogen temperatures under vacuum. The 
detectors are tuned to search for electrons emitted from 
the "Ge — Se double beta decay with a 2.041 MeV 
summed electron energy, which is the total kinetic en- 
ergy available for this particular reaction. 


The second approach to measuring the half-life of a 
double beta decay reaction is based on the cumulative 
effects of double beta decay over geological time. This 
geochemical method depends on the mass spectrometric 
determination of the stable double beta decay product 
that has accumulated in geologically old minerals. The 
geochemical method does not directly determine the de- 
cay mode as is possible in the direct determination 
method. lt indicates only the total decay probability of 
two neutrino decay X, and no neutrino decay A, A, + 
X. However, for double beta decay accompanied by the 
emission of two neutrinos, the theoretically predicted de- 

cay probability A, is much smaller than that predicted in 
neutrinoless double beta decay X „by perhaps six orders 
of magnitude (Doi et al. 1985), although the actual value 
depends on the mass of the neutrino. 


The geochemical method provides an upper limit to 
the double beta decay half-life because the gas retention 
age invariably post-dates the time of mineralisation; the 
direct detection method yields lower limits because of 
the possibility of background events. The direct detection 


r - process 
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method has the potential advantage of being able to dis- 
tinguish between the two double beta decay modes. In 
the case of the two neutrino double beta decay mode, the 
energy deposited in the detectors will show a continuous 
energy distribution up to the total transition energy, in 
the same manner as single beta decay. However, if 
neutrino-less double beta decay occurs, the sum of the 
energies of the two emitted electrons will be the total 
transition energy, and the energy spectrum will be a 
single narrow peak. The energy spectrum therefore pro- 
vides a clear, unambiguous signature of the nature of 
double beta decay. The high resolution, double focusing 
mass spectrometer at the University of Manitoba has 
been used to determine double beta decay energy values 
for WTe and “Te of 867.2 + 1.0 and 2528.8 + 1.3 keV 
respectively (Dyck et al. 1990). The total kinetic energy 
available for any double beta decay system can be ob- 
tained from the atomic mass table of Wapstra & Audi 
(1985). 


Geochemical half-life determinations of 
double beta decay 


The geochemical method of measuring double beta 
decay depends on the accumulation of decay products 
from this rare nuclear process in a natural mineral of the 
parent element. The method is most likely to succeed if 
little of the daughter element was incorporated when the 
mineral of the parent element formed. The likelihood of 
this separation occurring in nature is enhanced if the geo- 
chemical nature of the double beta decay product is radi- 
cally different from that of the parent element. Success 
also depends on the mineral being geologically old and 
having survived intact over a long period of time so that 
detectable quantities of the daughter have accumulated 
in the mineral. Hard, high temperature minerals have a 
favoured chance of surviving temperature/pressure 
variations to which the minerals may have been exposed 
during their geological history. 


Geochemical measurements of double beta decay 
have, until recently, been confined to the measurement 
of gaseous daughter products (by gas source mass spec- 
trometry), which have accumulated in old tellurium and 
selenium minerals from the following reactions; 


ae sn 
wie — BB = Me Zn 
Ha SIS a 2 


The double beta decay of 'Te to Xe is a favourable 
experimental situation because '*Te is the most abun- 
dant nuclide of Te, old telluride minerals of high concen- 
tration are available, and the mass spectrometry of Xe 
can be carried out with high sensitivity. Inghram & 
Reynolds (1949) first estimated the half-life of '"Te for 
double beta decay to be 1.4 x 10”y. Since that time, 
many "Te decay experiments have been carried out and 
a partial list of the results are given in Table 2. The 
current best estimate of the half-life of “Te is 8 10%y 
(Manuel 1991). Figure 2 (taken from Kirsten et al. 1968) 
shows the dominance of '*Xe in the mass spectrum of Xe 
extracted from an old telluride ore sample. Similar mass 
spectrometric determinations have been made for the 
double beta decay half-life of 'Te and “Se. Manuel 
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(1991) gives best estimates of these half-lives as 2 10%y 
and 1 10%y respectively. The results from these isotopic 
systems support the lepton-conserving, two neutrino 
mode of decay. 

Table 2 


Values reported for the half-life of "Te 


Reference Half-life 
Inghram & Reynolds (1949) 1.4 102 
Takaoka & Ogata, (1966) 324102 
Srinivasan et al. (1972) 25 10% 
Hennecke et al. (1975) 9:7 10% 
Kirsten (1983) 2.6 102 
Richardson et al. (1986) <1 102 
Kirsten et al. (1986) MONOS 
Manuel (1986) 22102 
Chiou & Manuel (1988) ZI 
Current best estimate 8 10% 


lon Current 


e je. 


128 130 132 134 136 


Mass Number A 


Figure 2. Isotopic composition of xenon extracted from tellu- 
rium ore (Kirsten et al. 1968). The horizontal lines indicate the 
maximum contribution of atmospheric xenon. 


The major analytical shortcoming of the Te/Xe and 
Se/Kr decay systems is that their daughters are gaseous 
and therefore may have suffered some leakage from the 
ore sample over the long time period since the minerals 
were formed. Minerals of tellurium and selenium are 
particularly susceptible to gas loss because they are soft, 
low temperature minerals that may have recrystallised 
and lost radiogenic Xe or Kr during any deformation 
event or high temperature phase subsequent to 
mineralisation. Richardson et al. (1986) have shown that 
gas retention ages are significantly less than 
mineralisation ages. Lin et al. (1986) discuss the prob- 
lems involved with gas retention in the mineral kitkaite. 
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Although there may well be uncertainties in the abso- 
lute values of '*Te, '®Teand ‘Se, there are compelling 
reasons why the ratio of the half-lives of '*Te, "Te for 
double beta decay should be determined. The great ad- 
vantage of the Te/Xe double beta decay system is that 
both the daughter products are noble gases and the ratio 
%Xe/'Xe, which is directly proportional to the ratio of 
the half lives of '*Te/™Te, can be determined with rela- 
tively high precision by noble gas mass spectrometry, 
thus eliminating many sources of systematic error. 


Pontecorvo (1968) has shown that the ratio of the de- 
cay rates of pairs of similar nuclei can be used to test 
various theoretical predictions because the ratio of their 
respective matrix elements should be near unity. In the 
case of Dirac decay, the ratio p, = A/M, is much 
smaller than the Majorana decay ratio p,, = 
where X is the probability of double beta decay. Thus 
even a small Majorana decay contribution would cause a 
drastic increase in the measured ratio p, = 'A,/™A, 
above the base level of Dirac decay p,,, (where È repre- 
sents the sum of the double beta decay mode compo- 
nents). 


The measurement of the half-life of the double beta- 
decay of '#Te is far more difficult than for ‘Te because 
the probability of double beta decay is so much smaller. 
However, a recent measurement of the double beta decay 
half lives of “Te and Te (Bernatowicz et al, 1993), has 
determined the daughter product ratio '*Xe/'WXe with 
greater confidence than previous measurements listed in 
Table 3 (Manuel 1991). The mean value for the ratio of 
half-lives in Table 3 is 4 10%, which compares favourably 
with the Bernatowicz et al (1993) value of 3.52 (+0.11) 107. 
The data of Bernatowicz et al (1993) give a limit to the 
effective Majorana mass of the neutrino of 1.5eV. 


Table 3 


Values reported for the ratio of the double-beta decay rate of 
Te relative to that of Te (Manuel 1991) 


Reference Dd See! as oe 
Hennecke et al. (1975) 6.3 107 
Kirsten (1983) FoTo 
Kirsten et al. (1986) 1.0 10° 
Manuel (1986) SA IO 
Takaoka & Sagawa (1988) <6 107 
Lin et al. (1988) 3.9 107 
Takaoka & Sagawa (1990) 322105 
Lee et al. (1990) ASADOS 
Current best estimate 4 107 


It took 20 years from the time that Pauli first sug- 
gested the existence of a neutrino to when Inghram & 
Reynolds (1949) undertook a determination of the half- 
life of double beta decay. However, a further 44 years 
were to elapse before the half-life of double beta decay 
was measured by solid source mass spectrometry. lm- 
provements in thermal ionisation mass spectrometry 
have occurred to such an extent over the past decade that 
it is now possible to conceive of the measurement of the 
half-life of double beta decay schemes which have a non- 
gaseous daughter product. 


In fact, Kawashima et al. (1993) have recently mea- 
sured the half-life of the double beta decay of *Zr from 
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the excess amount of “Mo found in old zircon samples. 
A solid source mass spectrometer was adapted to im- 
prove the sensitivity for Mo isotopic analyses by the use 
of ion enhancing media. The resulting half-life for the 
double beta decay of *Zr of 3.9 (+ 0.9) 10" y is lower than 
the measured half-lives of Te and “Se. Zircons have a 
particularly tight lattice structure and this property has 
been successfully exploited in U/Pb geochronology be- 
cause the radiogenic Pb has remained in-situ over exten- 
sive geological periods. Unfortunately the presence of U 
in the zircon samples generate spontaneous fission prod- 
uct Mo which make the measurement of double beta de- 
cay “Mo difficult to identify. Kawashima et al. (1993) 
were forced to assume that the spontaneous fission yields 
for U (which so far have not been measured) are iden- 
tical to those of the thermal neutron fission of “Pu to 
make the necessary corrections. The authors also report 
a “strange excess” in “Mo which they were unable to 
explain. However, despite these difficulties the work of 
Kawashima et al. (1993) has shown that isotopic systems, 
other than those with gaseous daughters, can be used to 
measure the half-life for double beta decay. 


An examination has been made of the various double 
beta decay schemes listed by Doi et al. (1985) which have 
non-gaseous daughters. The best candidate seems to be 
the decay of '212’5n >1212Te, This decay scheme has the 
advantage that two isotopes can experience double beta 
decay, old cassiterite samples of known age are avail- 
able, and the mass spectrometry of tellurium is well es- 
tablished (Loss ef al. 1990). Furthermore, this system has 
the great advantage with respect to the Zr/Mo double 
beta decay system in that cassiterite does not contain a 
significant amount of uranium, and in any event tellu- 
rium is close to the mass valley of nuclear fission where 
the cumulative fission yields are very low. In addition, 
cassiterite samples can be converted to the metallic form 
(McNaughton & Rosman 1991), and ion exchange chemi- 
cal techniques have been established to separate tellu- 
rium from other elements (Loss et al. 1990). Doi et al. 
(1985) predict a half-life for '*Sn which is approximately 
a factor of ten longer than that predicted for *Zr double 
beta decay. 


It is possible that accelerator mass spectrometry could 
be applied to the problem of measuring the half-life of 
double beta decay systems because of this technique’s 
ability to remove molecular interferences and the separa- 
tion of isobars by rate of energy loss or range measure- 
ments. This enables high sensitivity mass spectrometric 
measurements to be made as an alternative to conven- 
tional solid source mass spectrometry, although the ne- 
cessity for precise determinations of the isotopic compo- 
sition of the element concerned may present a problem. 


Conclusions 


It seems surprising that 65 years after Pauli’s postu- 
late of the existence of the neutrino, we still do not know 
some of its basic properties. Neutrinos were produced in 
huge numbers during the very early stages of the Big 
Bang, but because of their extremely low interaction with 
other particles they seem destined to inhabit the Uni- 
verse as lonely “ghosts of creation”. They are also pro- 
duced copiously from every Main Sequence star, and 
play an essential role in supernovae explosions (Clayton 
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1989). Some mass spectrometric evidence to support this 
latter phenomenon has recently been published (Maas et 
al. 1996). It is a strange quirk of Nature that double beta 
decay, which is an extremely rare nuclear phenomenon, 
may be the mechanism which will enable the properties 
of the neutrino to be better understood. The basic insta- 
bility of certain r-process nuclides contrasts with their 
present abundance distribution, which can only be ex- 
plained if double beta decay rates have an extremely 
small probability. 


The neutrino is the only fermion for which we do not 
know if it is different from its antiparticle and is there- 
fore a Dirac particle, or if it is identical with its antipar- 
ticle and therefore a Majorana particle. Double beta de- 
cay offers a novel means of testing Grand Unified Theo- 
ries of particle physics. In double beta decay two modes 
are possible - one in which two antineutrinos accom- 
pany the emission of two beta particles (the Dirac neu- 
trino) and the second in which no neutrinos are emitted. 
The clear signature of neutrinoless double beta decay is 
the emission of two electrons whose combined kinetic 
energy equals the mass difference between the parent 
and daughter nuclides. If the latter mode can be experi- 
mentally established, it may yield limits on the electron 
neutrino mass, indicate the degree of lepton conservation 
in elementary particle interactions and provide fresh evi- 
dence for neutrino oscillations. The large difference in 
half-life between the two modes of double beta decay 
offers a sensitive indicator of estimating the contribution 
of neutrinoless decay to the well established two neu- 
trino decay mode. 


A considerable body of experimental evidence for the 
half-life of double beta decay has now been established 
by both direct counting techniques and geochemical mea- 
surements, although the latter technique probably pro- 
vides an upper limit because of the high probability of 
gas loss from old tellurium and selenium-rich minerals. 
However, new geochemical methods based on the detec- 
tion of solid daughter products by sensitive solid source 
mass spectrometric techniques should resolve the gas 
loss problem. Existing evidence suggests that 
neutrinoless double beta decay does not occur to any 
significant extent, that the electron neutrino mass may be 
as small as 1.5eV and that lepton conservation is main- 
tained. However the evidence is fragmentary and 
plagued by experimental difficulties so that further work, 
especially on double beta decay, is required. 
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